Draft version August 24, 201 1 

Preprint typeset using I^TgX style emulateapj v. 1 1/10/09 



TESTING PHOTOMETRIC DIAGNOSTICS FOR THE DYNAMICAL STATE AND POSSIBLE IMBH PRESENCE IN 

GLOBULAR CLUSTERS 

Eva Noyola 

Max-Planck-Institut for Extraterrestrial Physics, Giessenbachstrasse 85748, Garching, Germany and 
Universitats-Sterwarte Miinchen, Scheinerstrasse 1, 81679, Munich, Germany 



o 

(N 

< 
(N 

< 

o 

6 



HOLGER BAUMGARDT 
Argelander Institute for Astronomy, University of Bonn, Auf dem Hiigel 71, 53121, Bonn, Germany and 
School of Mathematics and Physics, The University of Queensland, Brisbane, QLD 4072, Australia 

Draft version August 24, 2011 

ABSTRACT 

Surface photometry is a necessary tool to establish the dynamical state of stars clusters. We produce realistic 
HST-like images from N-body models of star clusters with and without central intermediate-mass black holes 
(IMBHs) in order to measure their surface brightness profiles. The models contain ~600,000 individual stars, 
black holes of various masses between 0% to 2% of the total mass, and are evolved for a Hubble time. We 
measure surface brightness and star count profiles for every constructed image in order to test the effect of 
intermediate mass black holes on the central logarithmic slope, the core radius, and the half-light radius. We 
use these quantities to test diagnostic tools for the presence of central black holes using photometry. We find 
that the the only models that show central shallow cusps with logarithmic slopes between -0.1 and -0.4 are 
those containing central black holes. Thus, the central logarithmic slope seems to be a good way to choose 
clusters suspect of containing intermediate-mass black holes. Clusters with steep central cusps can definitely 
be ruled out to host an IMBH. The measured rc/ri, ratio has similar values for clusters that have not undergone 
core-collapse, and those containing a central black hole. We notice that observed Galactic globular clusters 
have a larger span of values for central slope and rc/ r/, than our modeled clusters, and suggest possible reasons 
that could account for this and contribute to improve future models. 
Subject headings: globular clusters: general 
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1. INTRODUCTION 

Surface photometry has often been the initial tool to estab- 
lish the dynamical state of globular clusters. The fact that the 
observed radial density of most clusters appears to be well 
described by King models (King 1966) has been taken as evi- 
dence that these clusters are relaxed systems and that their dy- 
namical evolution is dominated by two-body relaxation pro- 
cesses. A natural consequence of two-body relaxation is the 
onset of core collapse, where the central density of a star clus- 
ter increasevwMeJhe core radius decreases (see section 1 .2 
of iNovola & Gebhardj (12006') for a detailed description of the 
process). Some clusters have been identified as having un- 
dergone core-collapse. These are cases with very concen- 
trated surface density profiles, showing steep central cusps 
with a central projected logarithmic slope of ^ -0.7 dCohnl 
[1980^, that depart from King-type cores. About 20% of the 
Galactic glo bular cluster population falls into this category 
(lTrageretalJll995h . 

Kinematical evidence for core collapse accompanied by 
tailore d models have been presente d for three cl usters: 
M15 (fDuiletaD '1997'), NGC 6397 dDrukieri IT995h . and 
M71 (Drukier et al. 1992). The expected velocity cusp has 
only been resolved and modeled for M15 (Baumgardt et alj 
UOOSb; McNamara et al. 2004; van den Bosch et al. 200C 
NGC 6752 has been considered to be a post core-collapse 
cluster by many authors, but different datasets and analysis 
methods find that it has a s r nall flat central core ( Lugger et al] 
[1991 iFerraroetalJlIOOl INovola & Gebh^ 120061) . The 
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size of the core might be consistent with models for gravother- 
mal oscillations ( Vesperini & Chernoff 199 4j). This cl uster 
shows a steep central velocity cusp (Dr ukier et al.ll2003h . but 
no tailored core-collapse model has been created for it. 

There are a variety of heating mechanisms that can drive 
energy into the core of a star cluster, causing it to expand, 
and thus preventing core-collapse. The effect of binary heat- 
ing by primordial binaries is the best studied mechanism to 
date (Gao et al. 1991; Vesperini & Chernoff 1994), although 
it has been proposed that most Galactic globular clust ers are 
not y et in the binary-burning phase of evolution (.Fregeaul 
120081) . The presence of stellar-mass black holes acting as 
an energy source, has recently been invoked to explain the 
distribut ion of core sizes in LMC and SMC globular clusters 
(iMackey et al. 2008). Mass loss by stellar winds during early 
time s of the cluster evoluti on also contributes to cluster expan- 
sion (iBaumgardt & Kroup a 2007; Hurlev 2007t iBastian et af] 
120081) . so clusters might expand considerably even if they are 
born with concentrated configurations. A recently proposed 
mechanism is velocity kicks imparted during white dwarf for- 
mation, which would also act as a heating me chanism for 
clusters with velocit y dispersions of a few km/s (iDavis et alJ 
lIOOllRegeau et alj|2009) . 

The presence of a central intermediate mass black hole 
(IMBH) of 100-10,000 Mq is another mechanism that can af- 
fect th e dynamical evolution of star clusters. Bahcall & WoU 
(1 19761) calculated the shape of the radial density profile for a 
single mass star cluster around a massive black hole. They 
predicted the formation of a steep central cusp with a log- 
arithmic slope of -1.75. iBaumgardt et al.l (l2004allbl) con- 
firmed these results based on direct N-body simulations. They 
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also showed that muhi-mass clusters with IMBHs are mass- 
segregated in their centers and that main-sequence stars have 
cusps that are significantly flatter than -1.75. They found that 
the IMBHs appear to produce shallow central cusps on the 
projected density profiles of bright main sequence stars for 
these clusters, with slopes of ^ -0.2, as opposed to steep 
power -laws (Baumgardt et al. 200 5|). iNovola & GebhardB 
(120061) and .Novola & Gebhardt (.2007h (hereafter caUed re- 
spectively NG06 and NG07) obtained surface brightness pro- 
files from HST images for Galactic, LMC, SMC and Fornax 
dwarf galaxy globular clusters. They found that about 20% 
of the globular clusters in their sample show central slopes in 
this intermediate range. 

The surface density profile shape can also be affected 
in the size of its core when a central black hole is 
present lTrenti et all (120071) estimated the value of r^/r/, (ra- 
tio of core radius to half-light radius) for N-body simulated 
star clusters containing central black holes. They used a den- 
sity averaged radius as a measure for the core radius. They 
found that the ratio tends to reach values around 0.3 for 
these cases, while the value is considerably smaller (< 0.1) 
for clusters without black holes. On the other hand, Hurleyl 
(|2007) finds similarly large rc/n, values for N-body simu- 
lations evolved including 0.5- 10% primordial binaries, but 
without a central black ho le. In this case, the Casertano & Hut 
method dCasertano & Hut 1985) was used to obtain the three- 
dimensional core radius. The different way in which the core 
radius was measured from the N-b ody simulations differs be- 
tween the two results. Recently, IVesperini & TrentH (120101) 
analyzed direct N-body models with and without IMBHs. 
They found that shallow cusps with logarithmic slopes as 
steep as -0.3 are present in various models, not only the ones 
containing black holes. The apparent discrepancies between 
different models using different analysis techniques stresses 
the importance of performing meaningful measurements on 
models so they can be properly compared with observational 
data. 

Direct dynamical evidence for the existence of central black 
holes using velocity dispersion measurements has been put 
forward f or three nearby glob ular clusters. M15 was the 
first case (iGerssen et al.l l2002l 12003). but alternative mod- 
els witho ut black holes were also shown to be good fits to 
the data (iBaumgardt et all l2003bl) . The latest detailed dy- 
namical measurement and model finds non-conclusive evi- 
denc e for the presence of a c entral black hole in this clus- 
ter (Ivan den Bosch et 10120061) . Gl, a large globular cluster 
in Andromeda, has stronger observational evidence to sup- 
port the presence of a cen tral black hole, frorn integr ated 
kinematical measure ments ([Ge bhardt et al.' '2003J, l2005h . as 
well as from X-ray (iPoolev & Rappaport 2006) and radio 
dUlvestad et al.l 12007 1) observations, but alternative scenar- 
ios hav e also been presented for this case ( Baumgardt et al. 
12003d) . Omega Centauri is the most recent case for which 
line-of-sight velocity dispersion measurements appear to sup- 
port the existence of a central black hole of 40,000 M© 
(jNovola et al."2008'), but proper motion measurements from 
HST images find differe nt results ( Anderson & van der Marell 
l2010t Ivan der Marel & Andersonll2oToir Evidence has also 
surfaced for intermediate mass black holes in extra-galactic 
disk galaxies based on X-ray observations. 

Ultra luminous X-ray sources (ULXs) have X-ray luminosi- 
ties higher than the Eddington limit for a stellar mass black 
hole. One of the possible explanations for this emission is 
that it comes from accretion onto an IMBH. For example. 
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the galaxy M82 contains an ULX source which is believed to 
host an IMBH based on the absolute brightness of the source 
(Matsumoto & Tsuru"1999^, 'Matsumoto et al.''2001), and its 
radio variability ( Strohmayer & Mushotzky 2003). The posi- 
tion of the X-ray source ap pears to coincide with the young 
dense star cluster MGG-11 dMcCradv et alj|2003i) . There is 
also the controversial case of the globular cluster RZ 2109 
in NGC 4472, which shows the first clear evidence for a 
star cluster hosting a black hole (Maccarone et al. 2007 ), but 
the si ze of the black hole is still under debate (Ze pf et alJ 
l2008h . One more interesting object is the X-ray source 
CXOJ03383 1.8-352604, associated with a globular cluster in 
the Fornax elliptical galaxy NGC 1399. Ilrwin et alj (120101) 
suggest the emission might come from a tidally disrupted 
white dwarf around an IMBH. 

In this paper, we create synthetic HST-like images from N- 
body simulations with and without IMBHs. We measure their 
surface brightness profiles as we would with observations. 
We provide an analysis of the detailed shape of central den- 
sity profiles for these models that helps understand the central 
state of Galactic globular clusters. We describe the N-body 
models in section 2, the synthetic images in section 3, data 
analysis in section 4, and discussion in section 5. 

2. N-BODY MODELS 

It is often challenging to make a direct comparison between 
the results of N-body models and observations because it is 
hard to take into account the sources of uncertainty of obser- 
vations. It is complicated to discriminate how much of the in- 
formation from the models would be available to an observer 
if the simulated object was on the sky at a realistic distance. 
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With the goal of making more meaningful comparisons, we 
take the output of N-body models and create realistic syn- 
thetic images from them. 

The simulations used here followed the evolution of star 
clusters with and without central intermediate-mass black 
holes. All star clusters contained 131,072 (128K) stars ini- 
tially and were simulated with the N-body program NBODY4 
jAarsethl 119991) using the GRAPE-6 computers at Tokyo 
University. Ste llar evolution was followed using the fit- 
ting formula of iHurley et al.l (1200 0V assuming a metallic- 
ity of Z=0.001. For most models, the initial density profile 
was given by a King Wq = 7 configuration but we also in- 
clude two models that started from a King Wo = 5 configu- 
ration. The deta i led des cript ion of the runs can be f ound in 
iBaumgardt et al.1 (l2003ah and lBaumgardt etal.1 ( 12005b . called 
BM03 and BMH05 respectively. 

Models mlt and m2t are models without IMBHs. The data 
was taken from the N=128K star runs in BM03b, who as- 
sumed a neutron star retention fraction of 1 0% and a mass 
range between 0. 1 and 15 Mq, according to a lKroupal (1200 lb 
mass function. Model m3t is a new model made for this pa- 
per, starting with N=128K stars distributed according to a 
iKroupa (2001) IMF from 0.1 to 100 Mq, and with an asu- 
umed neutron star and stellar black hole retention fraction of 
10% (.Pfahl et al.l2002|[?h . The retention fractions are assumed 
to be the same for simplicity, since there is still considerable 
uncertainty about these numbers. This simulation contains 
stellar mass black holes with masses up to 25 Mq until the the 
stellar-mass black holes have kicked each other out in two and 
three-body interactions at about T=12Gyr. All models with- 
out black holes are at a galactocentric distance of 8.5 kpc. 
Models mblt to mb4t come from BMH05, in this case the 
neutron star retention fraction was 15%, and the stellar mass 
range went from 0. 1 to 30 Mq, assuming the same mass func- 
tion as for the non IMBH cases. The models contain IMBHs 
of masses 0.2%, 0.5%, 1.0%, and 2.0% of the total mass of 
the star cluster (Mtot)- If we extrapolate the scaling laws for 
super-massive black holes in galaxies to the mass regime of 
globular clusters, the case with M,=0.5% Mjot would follow 
the Magorrian relation (Magorrian et al. 1998). Stars passing 
close to the IMBH were assumed to be tidally disrupted. We 
use the .Kochanek (^1992) formula for the disruption radius. 
We also performed one additional simulation of a star cluster 
with an IMBH (called mb5t in Table 1). For this simulation 
we overlayed four snapshots of the mb4t cluster at T = 1 1 
Gyr and continued the simulation for 1 Gyr with = 508.000 
stars. Given the large number of stars, no stacking was neces- 
sary for this cluster. We use this model to test if there is any 
effect from the stacking of close snapshots. 

For the non IMBH models, we created different snapshots 
in order to investigate the core-collapse evolution. Snapshots 
were taken at 1.0, 4.0,7.0,9.0, 10.0, 11.0, 12.5, and 16.0 Gyrs 
for the mlt case; and 2.0, 6.0, and 8.0 Gyrs for the m2t case. 
For the mlt model, core-collapse occurs at 12.5 Gyrs, while 
this happens at T=21.3 Gyrs for the m2t model and at T=20.5 
for model m3t. For the models containing IMBHs we use 
snapshots at different evolutionary times, all between 1 1 and 
12 Gyr, except model mb4t, for which we have earlier snap- 
shots. Information extracted from the models include mass, 
position, V magnitude, and temperature of each star. The de- 
tails of the created models can be found in Table [T] 

Given the initial number of stars and their mass function, 
stellar evolution, tidal evaporation, and disruption of stars by 
the IMBH, the final mass of the models is between 15,OOOM0 




Fig. 1. — Synthetic images for various N-body simulations. Tlie top and 
middle rows show a cluster without a central IMBH at different evolution- 
ary times of 1.0,4.0,7.0,9.0,10.0,11.0,12.5, and 16.0 Gyr (models mltl.O- 
mltl6.0 of Table 1). The bottom row shows clusters containing IMBHs of 
0.2%, 0.5%, 1.0% and 2.0% Mtot (models mbltll.5,mb2tll.5,mb3tll.5, 
and mb4tll.5 of Table l).The evolution toward core-collapse for the non 
IMBH case is clearly visible on the images in the upper and middle rows. 

to 45,OOOM0, which is only about 1-10% the mass of a typi- 
cal Milky Way globular cluster. Since the analysis performed 
in this work requires a large signal in the images, we had to re- 
sort to stacking snapshots separated by short periods of time 
around a given age for every model. For the models with- 
out an IMBH, we stack ^10 snapshots separated by 15 Myrs, 
while for the models containing IMBHs we stack 5 snapshots 
separated by 5Myrs. The ultimate goal is to have the same 
number of stars in the central region for every model. Mod- 
els with IMBHs are not subject to any external tidal forces, 
while the models without IMBHs are placed on a circular orbit 
around a Galactic tidal field and therefore undergo a stronger 
mass loss. In the end, the total number of stars present in our 
original lists is always around ^ 500,000. The total mass 
for the stacked models is ^ 220,000 Modot for the non-BH 
models and ^ 180,000 Mgdot for models with central BHs. 
Variations between individual models are under 10%. 

The total number of stars for each N-body model is given 
in column 6 of Table [T] The number of stars included in the 
synthetic images is given in column 7. As explained in detail 
in Section [3] this constitutes only ~ 50% of the original list 
due to brightness and radius cuts. The modeled clusters are 
also more extended compared to the Galactic clusters. With 
the goal of making the modeled clusters look more like dense 
Milky Way clusters, we also scaled the clusters down in size. 
We do this by dividing their coordinates by a common factor, 
which we chose to be a ^ 8 for the models without IMBHs, 
and a factor of 4 for the models with IMBHs. This scales 
all clusters down to a similar half-light radius (7-10 pc for 
non-collapsed cases), which is similar to that measured for 
Galactic globular clusters. 

3. CREATING SYNTHETIC IMAGES 

Our main goal is to create realistic images from the N-body 
models in order to perform the same type of analysis that we 
do on HST observations. The quality and size of the images 
is chosen to match that of the PC chip in WFPC2 or the HRC 
channel in ACS. In this way, we can make a proper compari- 
son with observed clusters contained on NG06. 

The procedure to create images is like the one described in 
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Fig. 2. — Star count projected density profiles from the input lists, wliicli contain more stars than those detected in the synthetic images. For every case, the 
magnitude bins are: Wmag < 16 (solid line), 16 <V„,a5 < 20 (dotted line), and Wmag > 20 (dashed line). The top row shows a case without an IMBH at evolutionary 
times of 1.0, 7.0, 10.0 and 12.5 Gyrs (models mltl.O, mlt7.0, mltlO.O, and mltr2.5). The bottom row shows cases containing IMBHs of 0.2%, 0.5%, 1.0% and 
2.0% Mjoj (models mbltll.5, mb2tll.5, mb3tll.5, and mb4tll.5). The evolution towards core-collapse affects the bright and intermediate bins, but not the 
faintest one. The presence of an IMBH affects the central slope of brighter bins, and the core radius of every stellar group. 



detail in NG06 and NG07. We use DAOPHOT (Stetson 1987) 
to add stars from a list of positions and magnitudes onto a base 
image. With the goal of including realistic background noise, 
we use as a base a WFPC2 image of a sparse field with the few 
present stars cleanly subtracted. We modify the base image to 
have a larger number of pixels than the PC chip on WFPC2, 
and we locate the center of the cluster at the center of the base 
image. The utilized point spread function (PSF) is obtained 
from observed data and it does not include variations across 
the chip. 

Since the center of observed clusters is not known a priori, 
we made a blind test, in which the center of the models were 
given an arbitrary shift in the three spatial coordinates, and the 
new center was calculated using the octants method described 
in detail on NG06. We choose a guess center and a radius, we 
count the stars present in eight 'pie slices' segments defined 
by the chosen center and radius and we calculate the standard 
deviation of the eight numbers. Using the same radius, we 
move to a new guess center and repeat the procedure several 
times around the initial guess center In the end we have a map 
of center locations and a standard deviation value associated 
to each of them. We fit a smoothing spline to the resulting 
surface and find the location of the minimum defined by the 
grid of guess centers, which we take as the true center For 
this procedure, we used every star in the list, which implies 
using many more stars than the ones that would be available 
to an observer Our goal is to test the method for a complete 
dataset, not to test the observed accuracy of the measurement. 



since this has already been tested in NG06 and NG07. The 
centers were calculated for three projections on the x-y, x-z, 
and y-z planes. In the end we found that the method is able 
to recover the center with an accuracy of O.Olpc ('^O.OOSrc). 
The tests performed in NG06 yield an error for the observed 
center location that corresponds to ^0.05rc. In general, the 
effect of measuring a density profile using the wrong radius is 
not necessarily to change the central surface brightness slope, 
but instead, a drop in the central measurement point is cre- 
ated. Seeing such a drop is actually an indication of having 
the wrong center, ? use this fact as a test for correct center- 
ing in their work, for example. Despite that drop, the slope 
of the other points up to the core radius is normally the same 
as the one using the correct radius. This is clea rly seen com- 
paring the profiles for omega Centauri betwe en iNoyola et aP 
(120081) and lAnderson & van der Marell (120101) . Despite using 
very different centers, and getting density profiles with differ- 
ent shapes, the slope of the profile between 15" and the core 
radius is consistent in both cases. 

The next step for making star lists suitable to be turned into 
images is projecting the stellar coordinates into a 2D distri- 
bution on the sky. We need to assume a fiducial distance that 
will affect both the coordinates and the magnitude of each star. 
The chosen distance for all cases was 5 kpc, which is on the 
near end of the distribution of distances for Galactic clusters. 
We choose this distance since it is adequate for our goal of 
obtaining high signal to noise images. After performing the 
geometrical projection and applying the distance modulus to 
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Fig. 3. — Surface brightness and star count profiles for models with IMBHs of 0.5% and 2.0% Mjqt at evolution times of 1 1.0, 1 1.5, and 12.0 Gyrs. The solid 
line is the N-body profile, the dotted lines mark the Poisson error for the N-body profile, the dashed line is the star count profile, and the filled points are the 
measured photometric points. The vertical scale is arbitrary. Uncorrected star counts underestimate central surface densities by a factor of 2-5, while photometric 
measurements are better suited to determine density profile for crowded fields, they always lie within the errors of the N-body profile. 



the star's magnitudes, bolometric corrections are performed to 
obtain the V-band luminosity of each star The correction is 
done taking i nto account the star's temperature following the 
procedure by iHurlev et alj ( 120001) . We create synthetic im- 
ages using DAOPHOT, which has a fiducial zero photometric 
point of 25 magnitudes, therefore, we eliminate from our list 
all stars fainter than that. These faint stars constitute ^ 20% 
of the entire list. At 5 kpc distance, Ipc radius is equivalent to 
41.25 arcsec. Assuming a pixel scale of 0.1"per pixel, this is 
equivalent to ^412 pixels. Taking into account the extra scal- 
ing factor mentioned in section |2] the synthetic images (1000 
pixels on the side) contain stars inside a radius of ^10 pc for 
each simulated cluster. Since we are interested in the central 
structure of the clusters, we choose the image size to include 
approximately 10 core-radii, and we exclude stars outside this 
radius. The final images end up including ^ 50% of the total 
number of stars in the simulated clusters. 

The results for a subset of the models can be seen in Fig. 1. 
For the model without an IMBH, it is clear that the cluster 
achieves a very concentrated configuration as it evolves to- 
wards core-collapse. On the other hand, the clusters contain- 
ing central IMBHs are less dense and have more extended 
cores. Once we have the synthetic images, we proceed to an- 
alyze them in the same way as we do with observed data. 

We count the number of detected stars inside the average 
core radius for our models. The average detected stellar den- 
sity in this region is ^2 stars/arcsec^. For comparison, this 



is an order of magnitude l ower than the centra l density de- 
tected for NGC 6388 by (iLanzoni et all l2007h . From the 
iNoyola & GebhardB ((2000) compilation, we located two clus- 
ters at different heliocentric distaces that have a similar central 
densities to our models, NGC 5634 (30 Kpc) and NGC 6541 
(7.5 Kpc) (.Harris 1996 ). 

4. SURFACE DENSITY PROFILES 

We measure surface density profiles for every synthetic im- 
age following the prescription described in detail in NG06. 
Using various DAOPHOT routines, we find stars and then per- 
form PSF-fitting photometry on them. DAOPHOT allows for 
the inclusion of noise when adding synthetic stars, therefore, 
even when we utilize the same PSF used to create the images 
for our photometric measurements, the subtractions are not 
perfect and are comparable to those in observed data. We have 
tested our measurement methods thoroughly using simulated 
images in NG06 and NG07. We know that we can measure 
the input centers within ~1 "for concentrated clusters, there- 
fore, we directly use the known input center for every image 
when we measure density profiles. 

The density profiles are obtained in two different ways: 
from integrated light and using star counts. A detailed discus- 
sion of the pros and cons for each method can be found in sec- 
tion 2.3 of NG06. For the first method, we use the magnitudes 
of detected stars to identify the brightest 2-3%, and we then 
proceed to mask them by giving them a value that excludes 
them from the integrated light measurement. For most stars. 
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Fig. 4. — Photometric points for various snapsliots of models witli and without IMBHs. The solid line is a smooth profile from the data points, the dashed line 
is a single-mass King fit. The vertical lines mark the different measured radii: black (thickest) is break radius, blue marks radius from King fit, and red (thinnest) 
marks FWHM core radius. 



we assign a masking radius of 3 pixels, which only eliminates 
the central bright region of stars, not the halo. Occasionally, 
if very bright stars are present near the center of the cluster, 
we use a larger radius to mask those. The haloes of the stars 
do contribute to the total light, but by masking the central part 
of the PSF disk, one prevents the giant stars from dominating 
the measurements. Obviously, the "contamination" effect is 
stronger in the very central regions in cases where there are 
many giant stars (like in post core-collapse). In our models, 
the integrated light fol lows the input profile very closely even 
in these cases. Also, iLiitzgendorf et alj ( 1201 1*) perform de- 
tailed PSF contribution estimations for ACS imaging of NGC 
6388. They conclude that the contribution of bright stars af- 
ter masking the central part of the PSF is under 10% for bins 
containing 10 pixels or more. Our bins are always larger than 
that. 

The number of detected stars is roughly 10% of the input 
stars, although it is worth pointing out that about 70% of the 
input stars are fainter than 20th magnitude. These stars make 
an important contribution to background light, but they are 
only detected as individual sources with low efficiency. As 
expected, the detection efficiency is close to 100% for the 
brightest stars (V„,aj<16), while the percentage declines for 



fainter stars, particularly closer to the center where crowding 
problems are worse. Once we have masked the 3% brightest 
stars, we measure integrated light by calculating the number 
of counts per pixel in various annuli using the bi weight, a sta- 
tistically robust estimator (Beers et al. 1990). As discussed in 
detail in NG06, this appears to be the optimal way to extract 
a density profile for stars with mass at or around the turnoff 
mass for an evolved cluster The choice of the sizes for the 
annuli is a tradeoff between obtaining the highest spatial res- 
olution and obtaining the least noisy profile possible. 

The second method we use to measure density profiles is 
star counts. From a star list, we construct a star count profile 
in the same annuli where we measure integrated light. This is 
done by estimating the number of stars per unit area, where 
every star has the same weight. As mentioned above, it is 
well known that in crowded field photometry, fainter stars 
are detected with decreased efficiency. The exact complete- 
ness fraction for a given brightness at a given radius, depends 
on the specific shape of each profile. Given that the surface 
brightness profiles are dominated by the brightest stars, we 
measure star count profiles only for the stars brighter than a 
given magnitude for each cluster, since this is the only way to 
make a meaningful comparison between the two methods. In 
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order to obtain formally correct star count profiles from im- 
ages, one must calculate the correct completeness correction 
factor for each brightness group in each image, which is very 
time consuming and outside the scope of this work. Uncor- 
rected star counts have been used to measure density profiles 
for star clusters recently (e.g. Lanzoni et al., 2008), so we 
feel that it is relevant to compare to such profiles. For the 
models containing IMBHs, the brightness cutoff always cor- 
responds to 16 magnitudes (slightly fainter than the turn-off 
point, equivalent to stars with 0.8 Mq); for the non IMBH 
models, the limiting magnitude changes with evolution time 
and is brighter than 16 magnitudes for every case except the 
most evolved case at 16 Gyrs. We use these limiting mag- 
nitudes to calculate a star count profile from the original in- 
put list, as opposed to the detected list, and we call this the 
'N-body profile'. We compare our measured profiles against 
this N-body profile, which can be thought of as the 'true' 
profile of the cluster, since it comes straight out of the en- 
tire model dataset. The precise limiting magnitude for each 
model is taken as the one for which the surface brightness 
profile matches the N-body profile in the region outside the 
core radius. 

We notice that every simulated cluster, with and without 
IMBHs, shows mass segregation, as can be seen in Fig. |2] 
where we compare profiles obtained from the input list for 
various brightness groups. As expected, the profiles for the 
brightest stars are more concentrated than for the intermediate 
and faintest groups. The faintest group almost always shows 
flat central densities, except for the case containing a 2.0% 
Mtot IMBH. As explained above, the solid lines in this fig- 
ure are taken as our 'N-body' profile for every model. 

We compare the measured surface brightness and star count 
profiles with the N-body profiles. This is shown in Figure [3] 
where we present three profiles for models with IMBHs of 
0.5% and 2.0% Mtot at different evolution times of 11.0, 
11.5, and 12.0 Gyrs. The limiting magnitude for the N-body 
profile and the measured star count profile is always the same. 
We note that the uncorrected star count profiles always under- 
estimate the density for the central regions, including at and 
around the core radius, while all three profiles agree very well 
at large radii. The N-body profile is sometimes noisy at the 
center, which is expected due to the small numbers of bright 
stars in that region. We show the Poisson noise for the N-body 
profile. As can be seen, for every case, the integrated light 
profile follows the N-body profiles very well at r > 30", and 
is as smooth as the N-body profile inside the core. The shape 
of the surface brightness profile is clearly dominated by the 
brightest stars, but the masking of the bright stars combined 
with the background contribution from fainter stars helps to 
make it smooth. It practically always lies within the Poisson 
errors for the N-body profile. 

Once we have obtained the ph otometric points for e ach 
case, we use a smoothing spline (IWahba & WangI ri990h in 
order to obtain a smooth profile for further analysis. Since we 
want to measure half-light radii as well as fit King profiles, 
we need to cover the complete radial extent for the clusters. 
Given that both surface brightness and star counts agree very 
well with the N-body profile at large radii, we extend the mea- 
surements using the N-body profile to the complete radial ex- 
tent of each modeled cluster We decide to truncate the star 
counts at the 0.2 pc width annulus for which we no longer de- 
tect stars. There might be stars present at larger radius, but we 
know that they are very few. The lower density limit we use 
is lower than what one could measure for observed Galactic 
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Fig. 5. — Photometric points for various snapshots of models without 
IMBHs but containing stellar-mass black holes. The solid line is a smooth 
profile from the data points, the dashed line is a single-mass King fit. 

clusters, where the field population already would dominate 
the measurements. In the end, we fit a smooth spline to a 
combination of our measured photometric points for the ra- 
dial extent of our images, and of the N-body profile at larger 
radii. 

5. ANALYSIS 

As mentioned in Section[T] two types of photometrical mea- 
surements have been proposed as possible diagnostics for the 
presence of IMBHs in star clusters, the central slope of the 
density profile, and the r^/r/, ratio. In this section we explain 
how we obtain both quantities for our simulated clusters. 

The measurement of the half-light radius (r/,) is straight- 
forward once we have the complete smooth profiles. We in- 
tegrate the light profile to get the total luminosity and take 
the radius at which the profile contains half the amount of 
light. The measurement of the core radius is more compli- 
cated since there are different definitions and ways to measure 
it for observations and numerical modeling. In this work, we 
explore three different ways to measure core radii that are nor- 
mally used for obse rved clusters. The first is the one used by 
iTrager et all (Il995h and iHarrisI (Il996h . whose results are the 
sources for most studies of large samples of Galactic globu- 
lar clusters. These catalogs define the core radius as the half 
width half maximum of the radial density profile (rd, from 
now on). This definition makes the radius resolution depen- 
dent when the profiles are not flat towards the center, since 
the closer to the center we measure, the brighter the central 
luminosity value becomes. A second definition is the one t hat 
comes from fitting a single-mass King profile (King 1962) to 
the density profile and taking the value of the fit for the core 
radius (which we call r^i,). The third definition is the one used 
in NG06, called break radius, and defined as the radius of 
maximum curvature of the density profile (called ri,). It can 
be understood as the turnover radius. 

As can be seen in Figure |4] and Table |2l the agreement be- 
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tween the three radii is good for models with central slopes 
between 0.05 and -0.05 (i.e., those with flat central cores). As 
expected, rd, is smaller than the other two radii for models 
with central cusps. It can also be seen that the King fits agree 
very well with the observed profiles for models with flat cores, 
while for the rest, the agreement of the King fit is good only 
outside the core radius, but the values of rck and start to 
diverge. Notice that the King fit does not describe the profile 
well towards the center. For models mltl2.5 and mltl6.0, 
the profiles are so steep, that we cannot measure a reliable 
minimum of curvature (r^). The King fit for these cases is 
a bad match for the entire radial extent, so, even though one 
can formally obtain a value for rck and rc/,, neither of them 
provide meaningful information about the density profile. For 
the cases in which we have three close snapshots, we notice 
that the deviation between the different radial measurements 
is of order 10% for ri,, 20% for r^-h, and 5% for r^k, but the de- 
viation between the three different types of core radii is larger 

The central surface brightness slope is obtained by calculat- 
ing the derivative of the smooth profile inside the core radius. 
This derivative is constant for r < r^. It is worth mentioning 
that the value is the same when we measure the slope of a 
linear fit to the photometric points in the same region. For 
the couple of very concentrated cases, models mltl2.5 and 
mltl6.0, where we cannot reliably measure a break radius, 
we take the central values of the derivative as the slope. 

If we try to use completeness uncorrected profiles instead 
of light profiles the value for r/,, ri, and rck does not change 
much, since it is the shape of the profile inside the core radius 
that changes, but not the turnover radius, rch on the other hand 
suffers a larger change since the value of the central density 
is lower. Obviously, the value of the central slope is very dif- 
ferent (flatter in general). If we try to construct a star count 
profile using only those stars that are detected with close to 
100% completeness, there are too few stars left and the pro- 
file becomes too noisy in the center to make any meaningful 
measurements. 

Fig. 5 shows the fit of single-mass King profiles to model 
m3t, which contains stellar-mass black holes. It can be seen 
that the fits are significantly more noisy in the center than 
those for models mlt and m2t, which is most likely due to 
the more stochastic heating of the cluster by a few black holes 
as compared to a core of neutron stars and white dwarfs. Ta- 
ble 2 shows that the derived photometric parameters are nev- 
ertheless still within the range seen for those of clusters with- 
out stellar-mass black holes, in particular the central surface 
brightness slopes are still all below -0.12. 

Once we have central surface brightness slopes and rc/ri, 
measurements for every modeled cluster we proceed to plot 
each point on a slope versus rc/ri, plane. We create a plane 
for each of the three measured radii. We find that using rck 
or rt gives very similar results, while using r^, does not give 
meaningful constraints for clusters without flat central light 
profiles, so we exclude this quantity from further analysis. In 
Figures 6 and 7 we show the location of our models on the 
slope versus / r/, and rck / fh planes respectively. The models 
span a range of central slopes from 0.18 to -1.00, but the only 
models that have slopes steeper than -0.5 are those that have 
achieved core collapse, while the steepest slope for a model 
containing an IMBH is -0.45. As shown in Figs 6 and 7, there 
are two models containing an IMBH that present a flat cen- 
tral core (model mb4t9.0 and mb4tll.5), otherwise, only the 
models containing intermediate-mass black holes show shal- 
low central cusps. Models without IMBHs show either a flat 
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central slope, or steep central cusps. Regarding the rb/ri, ra- 
tio, the cases that haven't reached core-collapse and started 
from a King model with Wo = 7 lie within a narrow range be- 
tween 0.15 and 0.35 and there is no clear distinction between 
these cases and those containing IMBHs in this respect. The 
cases that clearly separate towards large rc/rf, are those that 
started from King models with Wo = 5. The models contain- 
ing stellar-mass black holes lie close to the first group, but 
have larger r/,/r/, values. The two core-collapsed cases are 
placed at /r/, = 0, since we cannot formally measure a break 
radius for them. For the rck/ fh case, the actual values change, 
but the behavior is similar The only group of models that 
clearly separates from the rest in both plots are those with 
very steep central slopes and non-detectable turnover radius, 
which correspond to clusters that do not contain IMBHs and 
have undergone core-collapse. 

We overlay on both planes all the Galactic clusters in NG06, 
plus omega Centauri and Gl. For Gl we measure the cen- 
tral slope using the profile in Gebhardt et al. (2005), while rc 
and r/, values come from the analysis of Ma e t al.l (120071) . The 
first thing to notice is that the Galactic clusters occupy a larger 
area in the plane than the modeled ones. The two clusters for 
which there are kinematical indications of hosting an IMBH, 
Gl and omega Centauri have central density slopes shallower 
than -0. 1 , and their r^/ r/, values are different. Omega Centauri 
and Gl sit near the locus of our models, but both of them have 
more extreme values of rc/r/, than the models with IMBHs. 
Very concentrated clusters, like M15, which are assumed to 
have undergone core-collapse, do lie very close to the mod- 
els without IMBHs and long evolutionary times. It should be 
noticed that some individual Galactic clusters change location 
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Fig .6. — Central surface brightness slope versus the ratio of ri, / ri, . The full 
circles mark the location of models containing IMBHs, the full squares are 
for models without an IMBH, and the full pentagons are for models contain- 
ing stellar-mass black holes. Representative error bars for the central slopes 
(NG06) are shown on the top. The full triangles mark the location of 38 
Galactic globular clusters, while the open triangles are for 01 and omega 
Cen. Some individual globular clusters are labeled. 

from one plane to the other 

6. DISCUSSION AND CONCLUSIONS 

M15 was the first cluster for which the presence of a central 
massive black hole was kinematically investigated, mainly 
due to its concentrated central profile. It was only later when 
it became clear that a projected steep central cusp is not the 
expected behavior for a star cluster containing a black hole. 
This stresses the need to develop better diagnostics to dis- 
criminate suitable candidates for detailed kinematical mea- 
surements when looking for IMBHs. In this paper we have 
created realistic synthetic images from N-body models of star 
clusters with and without intermediate-mass black holes. We 
have analyzed these images in the same way we analyze HST 
data for a sample of Galactic globular clusters and we com- 
pare both datasets. We explore two quantities as possible di- 
agnostic tools for the presence of black holes: the central log- 
arithmic slope of the surface brightness profile, and the ratio 
of core radius to half light radius. We find that the r^/ r/, ratio 
cannot discriminate between models with and without black 
holes, as Hurley ( 2007) already found, but that the central log- 
arithmic slope can. N-body clusters without IMBHs show ei- 
ther flat central cores, or steep cusps if they have undergone 
core-collapse, while clusters containing IMBHs show shallow 
central cusp for all except two cases. 

We want to emphasize that when dealing with density pro- 
files of star clusters, saying 'core radius' alone is not enough, 
one has to specify how that radius is measured in order to 
compare its value to models or other observations. Histori- 
cally, the definition that we call rch is the most popular one, 
but we show here that this definition is only useful for clusters 
whose profiles have a flat central core. When the profiles have 
central slopes different from zero, rck or r/, appear to be more 



Fig. 7. — Central surface brightness slope versus the ratio of r^jf/ri,. The 
filled circles mai'k the location of models containing IMBHs, the filled 
squares are for the models without an IMBH, and the full pentagons are for 
models containing stellar-mass black hole. As in the previous figure, repre- 
sentative error bars for the central slopes (NG06) are shown on the top, the 
full triangles mark the location of 38 Galactic globular clusters, while the 
open triangles are for Gl and omega Cen. Some individual globular clusters 
are labeled 

suitable, although they can differ by up to a factor of two for 
the same cluster For density profiles with a flat central core 
and clear turnovers, all three definitions mark practically the 
same radius. 

There are various ways in which our N-body simulations 
are idealized compared to Galactic globular clusters. First, 
the number of stars and the central densities are lower than for 
real clusters. Increasing both quantities would increase the re- 
laxation time, which in turn would increase the evolutionary 
times. Including the presence of binaries could change not 
only the timescales but also the nature of the core contraction 
and expansion. Also, our analysis comes from images with a 
limited amount of signal to noise (a combination of number 
of stars and fiducial distance). It is likely that the compar- 
isons would be more meaningful using images with a larger 
number of stars. These issues have to be kept in mind when 
comparing to observed globular clusters. Despite the idealiza- 
tions in the models, the span of rc/r/, and central slope values 
seems to generally agree between our models and observed 
clusters. We note that the agreement between the two mod- 
els that have undergone core-collapse and the observed clus- 
ters that are suspect of having undergone the same process is 
very good. Having simulations with a larger number of stars 
would allow to analyze snapshots closer in time to fully ex- 
plore the process of core-collapse. There are some areas of 
Figs 6 and 7 containing observed clusters that our models do 
not populate. As mentioned in section [1] a variety of heat- 
ing mechanisms have been proposed for star clusters in recent 
years. Some of them, like mass loss or white dwarf kicks 
should affect most clusters; while others like tidal shocking 
or primordial binaries depend on the structure and evolution 
history of each cluster. A combination of including some of 
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these heating mechanisms and starting from a larger variety of 
configurations (a larger range of initial Wo) would likely pro- 
duce a better agreement between models and observations. 

Our results are in contrast with those of IVesperini & Tre"ntil 
(EOlO) since they find a number of models that present central 
shallow cusps without containing black holes. We think the 
reason for the difference between their result and ours lies on 
a combination of two things: on one hand, their models con- 
tain about 10% of the number of stars our models have. On 
the other hand, they count main sequence stars, which we find 
to be detected with a large degree of incompleteness in real- 
istic analysis, particularly in the center of rich clusters. Thus, 
we are tracing a different subset of stars when measuring den- 
sity profiles. We believe that the lower numbers of stars in 
their models produces noisier profiles that in turn can show 
shallow cusps due to fluctuations in the photometric points. 
This is illustrated by the fact that as soon as they use more 
particles (64K runs with combined snapshots), their central 
slopes before core-collapse times converge to shallower val- 
ues consistent with the ones we find for models without black 
holes. 

It is clear from figures 6 and 7 that a division can be made 
between clusters with and without black holes using only 
one of the two quantities that we have explored, the cen- 
tral logarithmic slope of the surface brightness profile. The 
rc/rh ratio cannot distinguish between cases with and with- 
out black holes, only between clusters that have undergone 
core-collapse and the rest. Clusters that have achieved core- 
collapse separate cleanly from the rest in both indicators, 
which leads to the exclusion of very concentrated clusters. 



like M15, as candidates for hosting an IMBH. None of our 
models reproduce central slopes between 0.5 and 0.65 and we 
observe a few Galactic globular clusters with those slopes. 
Since we are not able to follow the details of the evolution 
right before core-collapse due to the time intervals between 
snapshots, we cannot rule out the possibility that clusters in 
this stage could have intermediate slopes Even if clusters un- 
dergo such a phase, it is expected to be only for a very short 
time. Two of the 14 models containing IMBHs do not show 
a clear central shallow cusp. Even though it is impossible to 
draw statistical conclusions from such a small sample, we can 
say that the absence of shallow cusp does not imply the ab- 
sence of a central black hole. Therefore, some clusters with 
shallow cores might still be interesting candidates to follow up 
with kinematics. Finally, clusters with central slopes between 
-0. 1 and -0.45 are clear candidates for harboring central black 
holes since we can only reproduce shallow central slopes by 
including intermediate-mass black holes. We conclude that 
the central logarithmic surface brightness slope appears to be 
a good diagnostic tool for choosing star clusters candidates 
for harboring intermediate-mass black holes. 
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